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FUTURE STABILITY OF THE FLRW FLUID SOLUTIONS IN THE PRESENCE 
OF A POSITIVE COSMOLOGICAL CONSTANT 

TODD A. OLIYNYK 


Abstract. We introduce a new method for establishing the future non-linear stability of perturbations 
of FLRW solutions to the Einstein-Euler equations with a positive cosmological constant and a linear 
equation of state of the form p = Kp. The method is based on a conformal transformation of the Einstein- 
Euler equations that compactifies the time domain and can handle the equation of state parameter values 
0 < K < 1/3 in a uniform manner. It also determines the asymptotic behavior of the perturbed solutions 
in the far future. 


1. Introduction 

The Einstein-Euler equations for an isentropic perfect fluid with a linear equation of state and a 
positive cosmological constant are given by 

Qij = fij _ ( A > o), (1.1) 

Vif ij = 0 , ( 1 . 2 ) 

where G*- 7 is the Einstein tensor of the metric 

g = gijdx l dx\ (1.3) 

A is the cosmological constant, and 

T*- 7 = (p + p)v l v^ + pg 1 -* ( gijV l v■> = — 1) 

is the perfect fluid stress energy tensors with the pressure determined by the equation of state 

p = Kp (K > 0). (1.4) 

On spacetimes of the form [To, oo) x T 3 , the Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions to 
(11.1I) - (U.2I) represent a homogenous, fluid filled universe that is undergoing accelerated expansion. In [12], 
the future non-linear stability of these solutions under the condition 0 < K < 1/3 and the assumption of 
zero fluid vorticity was established by Rodnianski and Speck using techniques developed by Ringstrom 
in HU. Subsequently, it has been shown HUMS] that this future non-linear stability result remains true 
for fluids with non-zero vorticity and also for the equation of state parameter values K = 0 and I\ = 1/3, 
which correspond to dust and pure radiation, respectively. It is worth noting that the stability result 
for K = 1/3 established in [8] relies on Friedrich’s conformal method [5J [S] , which is completely different 
from the techniques used in mmm for the parameter values 0 < K < 1/3. 

In this article, we introduce a new method for establishing the future non-linear stability of FLRW 
solutions to inii-pi). The advantage of our method is threefold: (i) it leads to a relatively compact 
stability proof, (ii) the parameter values 0 < K < 1/3 can be handled in a uniform manner, and (iii) 
it relies on a conformal transformation that compactifies the time domain to a finite interval, which we 
expect will be useful for numerical simulating solutions globally to the future. A further advantage of 
our method is that it can be easily adapted to handle the case K = 0; we will present the details in a 
separate article. 

Our approach to solving the Einstein-Euler equations dD-ca shares much in common with the 
methods employed in hu innini in that it is based on a particular choice of wave gauge that brings 
the Einstein-Euler system into a form that is suitable for analyzing behaviour in the far future via 
energy estimates. We remark that, for us, a suitable form means a symmetric hyperbolic system of the 
type analyzed in Appendix [B] The key difference between our approach and HU SI HU is that we do 
not analyze the physical spacetime metric directly, but instead use a conformal metric as our primary 
gravitational variable. It is worth noting that our method does not rely on the conformal field equations 
of Friedrich Hi- 
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In order to describe our approach, we first must fix our notation. The spacetimes that we consider 
are of the form M = (0,1] x T 3 . We use ( x 7 ), (/ = 1,2,3), to denote (spatial) periodic coordinates on 
T 3 , and t = x° to denote a time coordinate on the interval (0,1] where the future lies in the direction 
of decreasing t. Unless stated otherwise, we use lower case Latin letters, i.e. i,j,k, to index spacetime 
coordinate indices that run from 0 to 3 while upper case Latin letters, i.e. /, J, A', will index spatial 
coordinate indices that run from 1 to 3. Partial derivatives with respect to the coordinates ( x l ) will be 
denoted by di = d/dx 1 . The fluid-four velocity v l is assumed to be future oriented, which is equivalent 
to the condition 

v° < 0. (1.5) 

As indicated above, we do not work with the spacetime metric m directly, but instead use the 
conformally transformed metric 

9ij = e ^9ij ■ ( 1 - 6 ) 

As has been demonstrated previously in G3 IS], a key technical advantage of the conformal approach is 
that it turns global existence problems into local ones via the compactification of spacetime. In our 
setting, this amounts to us being able to obtain solutions globally to the future by solving the conformal 
Einstein-Euler equations on a finite time interval. An additional benefit of the conformal approach is that 
it allows us to account for the leading order effects of accelerated expansion through a judicious choice 
of the conformal factor. 

Under the conformal transformation m, the Einstein equations m transform as 

G ij = T H e 4 *fij _ + 2(V' V J $ - - (2D$ + |V$| 2 g )g ij , 

or equivalently 


- 2 R Z1 = -4 ViV,$ + 4Vi$V,$ - 2 


□ <f> + 2|V$| 


1 - K 


p + A 


,2$ 


9ij - 2e 29 (l + K)pviVj, (1.7) 


where we employ the notation □ = g’^V,V, and |V<I>| 2 = g lJ V t d 1 Vjand we have introduced the 
conformal fluid four-velocity 


We note that the four-velocity future orientation condition (ED is equivalent to 


v 0 > 0 

where vo = gojV J ■ Unless otherwise specified, we raise and lower all coordinate tensor indices using the 
conformal metric gij and its inverse g' :l as appropriate. 

Letting T and T denote the Christoffel symbols of the metrics <jij and gij , respectively, the difference 
rf ? - — T* is readily calculated to be 

ft- - = 9 kl {9iN^ + ft/Vph - 

Using this, we can express the Euler equations ED as 

Vi fij = 67” j V yl> + g lm f lm g^^. (1.8) 

The wave gauge we employ is defined by the vanishing of the vector field 

Z k = X k + Y k (1.9) 

where 


X k := g ij = -d i9 ki + IgUgi&g*, 


and /i is to be specified. 

For the conformal factor, we choose 

while for /z, we choose 


Y k = -2V fc $ + fiSS, 

$ = - In (t) 

2A 


p = 


3 1 


With these choices the vector field ED becomes 

1 1 2 ( fcn , A . \ 


( 1 . 10 ) 
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Remark 1.1. Since we want the conformal factor e 2$ to become unbounded towards the future to account 
for accelerated expansion, the choice (11.101) for our conformal factor explains our choice of time orientation 
on the interval (0,1] and shows that future timelike infinity is located at t = 0. 

For use below, we introduce the background Minkowski metric 

7i ij = -^5i5l + 8\8 i J 5 IJ , 

the densitized three-metric 

fl /J = det (g K L)h IJ (1-12) 

where 

( 9u ) = 

and the variable 

q= 5 °°-r ? 00 -|ln(det( ff p O)). (1.13) 

Remark 1.2. While our choices for the coordinates and conformal factor are ultimately justified by the 
fact that they work, motivation for these choices can be seen by observing that 

is an exact solution to the vacuum conformal Einstein equations, i.e. El with p = 0, that is a member 
of the FLRW family of solutions and locally conformal to de Sitter space. 

With our notation and conventions fixed, we are now able to state our future non-linear stability result 
in the following theorem, which is the main result of this article. The proof is given in Section [3j 

Theorem 1.3. Suppose A > 0, e > 0, 0 < K < 1/3, k £ Z> 3 , g (/ £ H k+1 ( T 3 ), g l ^,po,Vj £ H k (T 3 ), 
Po{x) > 0 for all x £ T , and that the quadruple 

(g' 1 ,&t9' 3 ,P, vi) 1 1=1 = (9 %, g \ J , p 0 ,«?) (i-i4) 

satisfies the constraint equations 

(G l0 -T M )\ t=1 =0, and Z k | t= i = 0. (1.15) 

Then there exists a S > 0, such that if 

\Wo ~ r l tj \\H k + 1 + llfl , i 7 ||ff fe + llPollff* + IMIff* < 8, (1-16) 

then there exists a classical solution £ C 2 ((0,1] x T 3 ), p,vi £ C' 1 ((0,1] x T 3 ) to the conformal Einstein- 
Euler equations given by m and (El that satisfies the initial conditions (11.141) . the gauge condition 
Z k = 0 in (0,1] x T 3 , the regularity conditions 

g ij £ C°((0,1], H k+1 (T 3 )) D C°([0,1], H k (T 3 )) n ^((0,1], FT fc (T 3 ))nC 1 ([0,1], H fc - 1 (T 3 )) 


and 


p,q£C o ((0,l],F t (T 3 ))flC o ([0,l],F t - 1 (T 3 )), 


and the bounds 


A nn , , 2A 3 . , 6 L„ T j, , 3 ,f t 

q < ~9 (t,x) < —, 2^ < -g 00 (t,x) < -, -6 J < g IJ (t,x ) < J , 


and 


— < v 0 {t, x) < 



for all (t , x) £ (0,1] x T 3 . 
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Moreover, there exist a, £ H k 1 (T 3 ) with a(x ) > 0 for all x £ T 3 , such that the solution satisfies 


and 


n°ii 





\\g 0j {t)- (.v 0j + t', j )\\ H «-i 

< 

— t 2 ln(t), 

(1.17) 

II d t g 07 

(*)- 

- 2 f 


< 

t, 

(1.18) 

t~\g° 3 {t) 

-v 

oj )\\H^ + \\d I9 oj m H *-' 

< 

—t ln(t), 

(1.19) 


llq(f) 

— q(0)|U*, + ||^q(t)|| ff fe-i 

< 

t, 

(1.20) 

ll0 /J W 


< 

t, 

(1.21) 




II* 3{1+K) p(t) - cr\\ H k-i 

< 

2 ( 1-3 K) 

t + t u+o , 

(1.22) 




IM*) - w(o) 

< 

1-3 K 

tTT+ry 

(1.23) 


for all t e [0,1], where vi(0) = 0 if K ^ 1/3. 


We conclude with a few remarks: 


(i) For any specified 5 > 0, an open set of initial data satisfying the constraint equations (11.151) and 
the condition (11.161) can be constructed using a variation of the method employed in [Sj §3]. 

(ii) The FLRW family of solutions correspond to the solutions obtained from initial data (11.141) of the 
form 


{g l 0 J ,gl 3 ,P 0 , v°j) = (-Oo655g + b 0 S^jS IJ , —oiSqSq + hS^jS 1 - 7 , c, 0) 


where the constants a 0 , bo, c € K.>o and ai, 61 £ R are chosen so that the constraint equations (11.151) 
are satisfied. 

(iii) More detailed behaviour of the solution near t = 0 can be obtained by using the estimates (11.1711 - 
(11.231) from Theorem 11.31 together with another application of Theorem IB. II from Appendix [B] 

(iv) As discussed in Section 1.3 of [12], it should be possible using Ringstrom’s patching argument from 
m to generalize the stability result contained in Theorem II. 31 to other spatial topologies. 


2. The reduced conformal Einstein-Euler equations 

As discussed briefly in the introduction, our proof of Theorem 11.31 relies on expressing the conformal 
Einstein-Euler equations, given by (ED and ED, as a symmetric hyperbolic system of the type analyzed 
in Appendix [B] We achieve this through the use of a wave gauge defined by the vanishing of the vector 
field (11.91) in conjunction with a suitable transformation of the field variables. We present the details of 
this procedure in Sections 12.11 and 12.21 


2.1. The reduced conformal Einstein equations. Following H. Friedrich’s [4] generalization of the 
gauge reduction method of Y.Clioquet-Bruhat !2] Ch. VI, § 8 ] that, crucially, allows for gauge source 
functions, we replace the Einstein equations ED by the gauge reduced version: 


—2i^ ? 2 ^-\~Aij Zk — — AViW + 4Vi4 ) V ? 4> 


^ 2 


□ $ + 2|V$ 


1 - K 


p + A e 


,2$ 


gij — 2(1 + K)e I<p pv 


l U J 


where 

A*h = -X^ +Y^. 

We will refer to these equations as the reduced conformal Einstein equations. 

Since the gauge condition Z k = 0 propagates, we can, with out loss of generality, assume that it holds, 
or equivalently that 

X k = -Y k , 

which, in turn, implies that 

X^H = Y° =^g 00 +^j. 

A straightforward calculation using (11.101) , (11.111) , (12.11) and the identity 

v {i x j) =-\d t g ij C x j = s J 0 ) 


( 2 . 1 ) 
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then shows that reduced Einstein equations are equivalent to 

-2 R ij + - X'X* = |^ (g 00 + - ^g 0K 5^ 


2(1 + K)^v i v j , 
t Z 

( 2 . 2 ) 


which we note can be expressed as 
-g kl d k dig ij - Q ij (g,dg) = ^d t g^ - ^ (g 00 + ^<Sj5g - 

^ (<?°° + |)-(i- *0^ - 2(1 + tf)£i;V (2.3) 


where the Q lJ are quadratic in dg = (<9fcg !J ) and analytic in g = (g y ) on the region defined by det(g* 7 ) < 0. 

The first step in transforming (12.31) into the desired form involves parameterizing the three-metric g IJ 
using densitized spatial metric (11.121) and the determinant variable (11.131) . Using the identities 

-g U d k d l ( 9 IJ ) = det (g PQ )ic I L J M (-g lk d l d k g LM ) - g^d^detigpQ^Ci^dtg^ 


and 

d t Q IJ = det (g PQ )^C I L J M d t g LM , 

where 

rIJ _ si sJ 1 - IJ 

LM — d L°M ~ 3 gLMg 3 

we see from m that q ij and q satisfy 

-g kl d k dig IJ - Q IJ {g,dg) = ^9 t g IJ - 2(1 + K) det{g PQ )^C I L J M -^v L v M (2.4) 

and 

-g kl d k diq-Q(g,dg) = |^d t q-2^ g + 3 ^ + (l-K)(^-g 00 ^J-^ + 2(l + K)(^g IJ v I v J 

(2-5) 

respectively, where Q IJ and Q are quadratic in dg = (d k g u ) and analytic in g = (g u ) on the region 
defined by det (g IJ ) > 0 and goo < 0. 

We proceed by writing the wave equations (12.31) (i = 0), (12.41) and (12.51) in first order form using the 
following definitions: 


u 0 - 7 

gOj _ gOj 

21 

(2.6) 

r. 

0 0 

to. 

= a , 9 »3_3( 9 «-« > 

2 1 

(2.7) 

0? 

U/ 

= d l9 °\ 

(2.8) 

u IJ 

= 0 /J - s IJ , 

(2.9) 

I J 

u k 

■"i 

O 

II 

(2.10) 

u 

= q, 

(2.11) 

Ufc 

= d k q. 

(2.12) 


and 


In terms of theses first order variables, it is not difficult to verify that the wave equations (12.31) (i = 0), 
dUD , and (ESI are equivalent to the following systems of symmetric hyperbolic equations: 


A k 8 k 
A k d k ( 





+ F, 


2 g°° 
t 


n 



+ 


(2.13) 


(2.14) 
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and 


A k d k 



,oo / u o\ 

n tJ +f 


(2.15) 


where 
A° = 

1 K 


A K = - 


P = 


A = 


n = 


f = 


5 = 



+ 4u 00 « + u oi ) - 8u 00 u 0Z + 6u aI u°j l - (1 - K)£g 01 - 2(1 + K) 

0 
0 

'Q lm + 4u 00 Uq M - 2(1 + K) det (g PQ )?£v R v s 

0 

„00 ,,LM 

~9 u o 


4rV°V 1 ^ 


(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


and 


f = 


Q + 4u 00 u 0 - 8(u 00 ) 2 + (1 - K) 


g 00 ^ $ + 2(1 + K) 
0 

_ n 00 u LM 

y u 0 


gRSV R V S - V°V° ) & 


(2.23) 


It is important to note that the term -A that appears in (I2.21|) - (12.23|) is not a singular term. This is 
because we can bound the density variable £ defined by (12.251) in the following section. Expressing £ in 
terms of it is clear that this term is regular. 

Another important point is that the first order variables (I2.6M2.12I1 are completely equivalent for 
t > 0 to the metric g lJ and its first derivatives d k g 13 . Consequently, we lose nothing by parameterizing 
the gravitational field in this fashion, and gain by having transformed the reduced conformal equations 
into a symmetric hyperbolic system of the forn{3 considered in Appendix [B] 


2.2. The conformal Euler equations. For the conformal Euler equations m, we use the formulation 
from |10l §2.2]; see, in particular, equation (2.55) of that article. In this formulation, given the conformal 
factor (11. 101) and the linear ecmation of state the conformal Euler equations m are from the 

computations presented in [lOl §2-2] readily seen to be given by 

(£) = i Br («,) - H < 2 - 24 ' 


lr To be precise, we still need to make the trivial coordinate transformation t > — t to obtain the form considered in 
Appendix [B] 
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where 


and 


P = i 3(1+K) /O c e (1+K)C , p c S M>o, 


= Vo(g ,v L ) ■■= 


-g 0I vi - yj ( g 0I vj ) 2 - g 00 (g IJ v I vj + 1 ) 


9 


oo 


V 1 = v , (g kl ,v L ) := g Vj + g M Vo(g ,v L ), 
J 




IJ 


JO 


v 


9 n i 

V v 


B v° (o -=^V JI 1 ’ 


Vq 


7r = 


0 0 


0 5j ' 


L k j = 6 k j- 


Vo 


J o> 


V I V J . 9 00 

Mij = gu - goj - 9i o + 7—To ViVj, 

v 0 v 0 {vo) z 


B k = 


K %L° m V mj \ 

£V li L° l V li M lm V mj ) 

1 f Kv k 


KL«V MJ 


v o \ KV li L% V li M lm V mj v k 


1 


-KL k L T L kl v l - KL\^d k g l - 


jj ___ | J-/ kl ~ ~^L dg lr 

v° \-v LI M LM v k T%v l - V LI M LM v k -§^d k g w ‘ 


(2.25) 

(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 

(2.32) 

(2.33) 

(2.34) 

(2.35) 


3. Proof of Theorem 11.31 

We begin by fixing A > 0, k € Z> 3 , K e (0,1/3], e > 0 and R > 0. Assuming g 0 o < 0, g 00 < 0 and 
det(g IJ ) > 0 , we then observe that the relations 


pa°p < pap and iL4°n = i (-2 ff 00 n) 


follow directly from (12.161) . (12.181) and ( 12 . 201 ) . 
Next, setting 


< A 0 0 0 

A° = I 0 A 0 0 


f A 

0 

0 


0 0 A 0 

0 


0 


and 


a straightforward calculation shows that 


- 2 g 00 l 

0 

0 


-2 g° 

/P 

0 

°\ 

0 

n 

0 * 

\o 

0 

n / 

LM 
*0 ! 

U L J M 

,u LM 


while we see from (12.161) that 


[A°,P] = [A,P] = 0, 
A 0 = A°(t,tu,P‘ L u). 


Using m and we compute 

D u A 0 • (A°) _1 APu = i[D 2 A°(i,tu,P J -u) • (A°) _ 1 AP: 

Next, it is clear from (12.291) and (12.301) that 


:= I —: 


(3.1) 


(3.2) 

(3.3) 

(3.4) 


[B, 7 r] = 0, 


(3.5) 
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and that 
Setting 


nB = 0, 


if K = 1/3. 


(3.6) 


£ = (u,C ,vj), 


and evaluating B° and B at £ = 0, we find that 


B 0\ 


l£=o = 

From these expressions, it is clear that 
7rB°| 4=0 7r = 

We also observe from (12.3311 that 


o L) and Bk = o= 


I 0 

0 1-3 K) 


1 


1 - 3 K 


7tB| c =o = (0, if Kg (0,1/3). 


(3.7) 


(3.8) 


'ir L B°{t,£ } )ir L = K (tt 1 - := II —7r), 

which, in turn, implies after differentiating that 

TT 1 - [D^B Q {t,C) ■ t^TT 1 = 0 

for all (i, £) in the domain of definition of B° and all variations <5£. 

From the evolution equations (12.131) . (12.141) . (12.151) . and (12.241) . the relations (13.11) . (13. 21) . (13.41) . (13.51) . 
m , » and (13.81) . and Remark |B.21 it is clear after performing the trivial time-coordinate transfor¬ 
mation 1 1 —y —t that £ satisfies a symmetric hyperbolic system of the type to which we can apply Theorem 
IB.11 from Appendix [Bj Doing so guarantees the existence of a 6 > 0 such that if 

(3.9) 


(ff’ J ,%* J ,C,«j)| t=1 = ( 9 o ,9 t g l 0 J Xo,Vj)\ t=1 € H k+l ( T 3 ) x H*{ T 3 ) x ff*(T 3 ) x H k (t d ) 


jk /qr»3 \ 


jk /r[ri3 \ 


rk /nr<3\ 


and 


Pc S 


is chosen so that p c and the corresponding initial data £|t=i = £o> which can be computed via the formulas 
11.121) . (11.131) . (12.6[l - (12.12[l and (12.251) . are bounded by 

ll£o||ff fc + Pc < S, 


then there exists a map 

£ e C 1 (( 0 , 1 ] X T 3 ]) D C°((0,1],i7 fc (T 3 )) D C°([0,1], Ff fc_1 (T 3 )) 


such that 

IICIlL“([0,l],ff fc ) < R, 

£ is the unique classical solution to (12.131) . (I2.14|) . (12.151) and (12.241) on (0,1] x T 3 that agrees with 
the initial data at t = 1, i.e. £|t=i = £01 the metric g l] and conformal four-velocity component vo, as 
determined by the formulas (11.121) , (11.131) , (12.61) , (12.91) , (12.111) and (12.261) , satisfy 


A 


<-g w (t,x)< 


2A 3 . , 6 

M <-»»«, x)< p 


3 ./j 


-S IJ <g IJ (t,x)<^S 


and 


for all (t, x ) € (0,1] x T 3 . Moreover, for given e > 0, it follows, by choosing R > 0 small enough, that the 


components of £ satisfy: 

\\u°o j (t) - u^'(t) + 2u° J (0)||#fc-i < t, (3.10) 

ll u o J (t) +u°- 7 (t)|| /ffc -i + llu^lljjfe-i < — tln(t), (3.11) 

||u /J (f) - u /J (0)|| ff fc-i + ||uj/(t) - u"(0)|| ff fe-i + ||u"(t)|| ff ,-i < t, (3.12) 

||u(t) — u(0)|| ff fc-i + ||u/(f) — U/(0)11 jjk 1 + ||9tu(t)|| H fc-i < t, (3.13) 

IIC(t)-C(0)||^-i + (3.14) 










































FUTURE STABILITY OF THE FLRW FLUID SOLUTIONS 


9 


and 

1-3 K 

|| Vi(t) - w/(0)|| H fc-i < t o+u (3.15) 

for 0 < t < 1, where 

ti/(0) = 0 if K G (0,1/3). 

Using (13.10p ~ (l3. 151) . it follows directly from the formulas (11.121) . (11.131) . (I2.6D - (12.12I) . and (12.251) that the 
metric and fluid density satisfy: 

II g° 3 (t) - ( if 3 + t'y 3 )\\ H k-i < —t 2 ln(t), 

II dtg° 3 (t) - 2t~ 1 {g° 3 (t) - if 3 ) + 7 J || ff fc-i < t, 

II d t g° 3 {t) - t~ 1 {g° 3 {t) - rf 3 )\\ H *-i + \\dig° 3 \\ H k-i < -tln(t), 

l|q(*) - q(o)||/f* + ||9tq(t)||ffk-i < t, 

II 3 U (t) - 9 IJ W\\ H k + ||<9 t g /J (f)||// fc -i < t, 

and 

in , 2 ( 1 - 3 K) 

||t p{t) — a\\ H k-i <t + t < 1+e > , 

where a,') 3 G U fc-1 (T 3 ), and a(x) > 0 for all x G T 3 . 

Since £ satisfies (12.131) . (12.141) . (12.151) and (12.241) on (0,1] x I 3 , it is clear from the results of Section 
12.II that the triple [g 13 , p = £ 3 ( 1 +- K ’) p c e^ 1+K ^,vi) determined by ^ solves the reduced conformal Einstein- 
Euler equations given by m and (12.21) on (0,1] x T 3 . This fact together with the assumption that the 
initial data (13.91) satisfies the constraint equations (11.151) allows us to conclude via the Proposition on pg. 
540 of [2] that the triple ( g 13 , p = p c e^ 1+K ^,vj) also satisfies the full conformal Einstein-Euler 

equations, given by and m, and the gauge condition Z k = 0 on (0,1] x T 3 . This completes the 
proof of Theorem 11.31 
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Appendix A. Calculus inequalities 


In this appendix, we collect, for the convenience of the reader, a number of calculus inequalities that 
we employ throughout this article. The proof of the following inequalities are well known and may be 
found, for example, in the books 00 > 0 and m- Before stating theses calculus inequalities, we first fix 
our notation and introduced a few definitions that will be needed in this appendix and the following one: 

Coordinates, indexing and derivatives: Lower case Latin indices, e.g. i,j,k, will run from 0 to n, 
while upper case Latin indices, e.g. /, J, K, will run from 1 to n. We use x = (x 1 ) to denote the standard 
periodic coordinates on the n- Torus T” and x° = t to denote the time coordinate on intervals of the form 
[To,Ti) where Tq < T\ < 0. Furthermore, we use lower case Greek letters to denote multi-indices, e.g. 
a = (or, 02 ,..., a n ) G Z" 0 , and employ the standard notation D a = d^d^ 2 ■ ■ ■ 9“ n for spatial partial 
derivatives. 


Inner-products and matrix inequalities: We employ the notation 

(£IC) = £ T C, 

for the Euclidean inner-product on M. N , and we denote the L 2 inner-product on T” by 

(u|u) = f (u(x)\v(x)) d n x. 

J T" 

For matrices A,Bg Mat x at, we define 

A< B «=» (£|A£) < (£|.B£), V^gR^. 


Constants and inequalities: We use that standard notation 

a<6 
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for inequalities of the form 


a<Cb 


in situations where the precise value or dependence on other quantities of the constant C is not required. 
On the other hand, when the dependence of the constant on other inequalities needs to be specified, for 
example if the constant depends on the norm ||u||l°°, we use the notation 


C = c 


Constants of this type will always be non-negative, non-decreasing, continuous functions of their argu¬ 
ments. 

With the preliminaries out of the way, we now state the calculus inequalities: 

Theorem A.l. [Holder’s inequality] If 0 < p, q,r < oo satisfy l/p+ 1/q = 1 /r, then 

\\uv\W < IMMMk* 

for all u £ L p ( T") and v £ L q (T n ). 

Theorem A.2. [Sobolev’s inequality] Suppose 1 < p < oo and s £ Z >n / p . Then 

IM|l~ ;$ IMIw^p 


for all u £ W s ’ p ( TT). 


Theorem A.3. [Product and commutator estimates] 
(i) Suppose 1 < pi , p 2 , 91,92 < oo, s = |o| £ Z>i, and 



Pi P2 q i 92 


Ther Q 


1 

r 


\\D a (uv)\\ L r < \\D s u\\ L Pi\\v\\ L n + \\u\\ L P 2 \\D s v\\ L <u 


and 


|| D a (uv) - uD a v\\ L r < ||Du||lpi \\v\\ w *-i, qi + \\Du\\ W s-i, P 2 \\v\\ L n 
for all u,v e C°°(T”). 

(ii) If Si, S 2 , S 3 € Z>o, si, S 2 > S 3 , 1 < p < oo, and si + S 2 — S 3 > n/p, then 

||nu||w s 3.p < |M|vv»i’p|M|w s 2,p 
for all u e W S1 ’ P { T") and v e W S2 ’ P ( T"). 

Theorem A. 4. [Moser’s estimates] Suppose s £ Z>i, 1 < p < 00 , |a| < s, f £ C S (R), /( 0) = 0, and V 
is open and bounded in R. Then 

\\D a f{u) ||lp < C(\\f\\ Cs(V) )(l + ||n||r~)l|w|l^.p 

for all u £ C°(T”) (~l L°°(T") (~| W s ' p { T") with u(x) £ V for all x £ T". 

Lemma A.5. [Ehrling’s lemma] Suppose 1 < p < 00 , so,s,si £ Z>o, and sq < s < si. Then for any 
e > 0 there exists a constant C = C(e) such that 

IMIw s ’P < e||u||rv' s i’P + C'(e)||'u||w' s o.p 


for all u £ W Sl ’ p (T n ). 


2 Here, we are using the standard notation ||Zt s u||^ p = X]| a |=s 
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Appendix B. Symmetric hyperbolic systems 


In this appendix, we analyze the initial value problem for symmetric hyperbolic systems of the form: 

1 


B l d iU = -BPu + H + G in [T 0 , T ± ) x T\ 


(B.l) 

u = u 0 in {To} x TP, (B.2) 

where 

(i) To < T\ < 0, 

(ii) IP is a constant, symmetric projection operator, i.e. P 2 = P, P T = P and dj? = 0, 

(iii) v = v{t,x) is a K M -valued map, u = u(t,x), G(t,x ) and H = H(t, x,v,u) are R^-valued maps, 
H £ C°([T 0 ,0], C'°°(T" x M m x R w )) and satisfies H(t, x,v,0) = 0, 

(iv) B 1 = B‘(t, x, v , u), and B = B(t , x, v, u ) are Mjv x jv-valued maps, and B 1 , B £ C°([To, 0], C'°°(T" x 
R m x 1^)), B° £ C' 1 ([T 0 ,0],C' oo (T" x R m x M N )) and they satisfy 

(B i ) T = B i and [P, 23] = 0, 

(v) there exists constants k, 71,72 >0 such that 


— 1 < B°(t,x,v,u) < —B(t,x,v,u) < 72 1 

7l K 


for all ( t,x,v,u ) £ [To, 0] x T” x K M x M , 

(vi) 

P ± B°(t, x , v, P^ ujP = P B°(t, x , v, P^wjP^ = 0 
for all (t,x,v,u) £ [T o ,0] x T” x R m x where 

P 2 - = I - P 

is the complementary projection operator0 

(vii) and there exists constants 0 , /3i, fc, /? 3 , P 4 > 0 and oj > 0 such that 


(B .3) 


(B.4) 


I_L [D u B°(t, x, v, u) ■ (B°(t,x,v,u)) 1 £?(f,:r,u,it)Pw]P J ‘| op < \t\9 + 
(P 2 - [ D u B°(t , x, v, u) ■ ( B°(t , x, v, u))~ 1 B{t , x, v, u)Pu]P| op < \t\8 + 


2pi 


[D U B°(t, x,v,u) ■ (B°(t,x,v,u)) 1 B(t, a;,i;,u)Pu]P J “| op < \t\6 + 


u> + IP^mI 2 

y/u + IP-Luj 2 
2/3 3 . 

\fui + |P^m | 2 


^| 2 , 


and 


s, [D u B () (t, x, v, u) ■ ( B°(t , x , v, u)) 1 B(t , x, v, w)Pu]P| op < \t\6 + /?4 


(B.5) 

(B.6) 

(B.7) 

(B.8) 


for all (t,x,v,u) £ [To,0] x T” x R M x Mr. 

Before proceeding, we make some observations. First, we note that the bound (IB.31) implies that B° 
and B are invertible and satisfy the matrix operator bounds 


|(T u (t,x, v,u)) | 0 p < 71 and \B(t,x,v,u) | 0 p < 


7i 


for all (t,x,v,u) £ [T 0 ,0] x T" x R M x M N . Second, a straightforward computation shows that the 
condition (IB. 41) is equivalent to 

P B°(t, x, v , P^uj-^P- 1 = P ± B°(t, x, v, P^u^P = 0 (B.9) 

for all ( t,x,u ) £ [Tq,0] x T" x 1^. We also note that it follows immediately from the bounds (IB. 31) that 
there exists constants R, 71,72 > 0, where 71 < 71 , k < k and 72 < 72 , such that 


TT-P < P B°(t, X , V, Zt)P < rrPBCt, X, V, u)P < 72 P 
7l K 


(B.10) 


for all (t,x,v,u) £ [Tq,0] x T” x M ai x 


3 In other words, P x .B 0 (f, x, v, it)P = P ± [B°(t, x,v,u) ■ Pu]P and PB°(t, x, v, njP 1 - = P [B°(t,x,v,u) ■ Pu]P for matrix- 
valued maps v, u ) and B°(t, x, v,u) with the same regularity as B° (t, x,v,u). 
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Theorem B.l. Suppose k £ Z >n / 2 +i, Mo £ H k (T n ), v,G £ C° ([To, 0], H k (T n )), assumptions (i)-(vii) 
are fulfilled, and the constants {/?, k, 71} satisfy 0 < ft/ 71 < k where ft = Xw=i 84 - TTien f/iere exists a 
T * s (To,0), and a unique classical solution u £ C 1 ([To,T*) x T") t/iaf satisfies u £ C°([To,T*), H k ) D 
C 1 ([To, T*), H k ~ 1 ), the energy estimate 

ll u (0llN + ll^llL“>([T 0 ,o),ff fc ) ~~ f “II^MNlN <Ce c<yt T °' > f||M(To)||j /fc + ||G||N([T 0 ,o),.ff fe )) 

Jtq 

for To <t <T* where 

C = C'(||M|| L c»([ r 0 )T< ) )H fc), ||M|| iO o([ T 0 | 0 ) iff fc), ||G||N([T 0 ,o),ir fc )! ||3tM|| L oo([ ToiO ) i iffc-i),0,7 1 ,72, k,ft,u)), 
and can be uniquely continued to a larger time interval [Tq,T*) for some T* £ (T*,0) provided that 
IMlL“([T 0 ,T.),wn.“) < OO. 

Moreover, for any R > 0, there exists a 

5 = 6(R, |M| L oo([ Toi 0 ) iH -fc), ||9 t M|| i oo([ To O ) jff fe-i),0,7i,72,/3,K,w) > 0, 

such that if || II ^r fc + IIGIU °°([T 0 ,o),H k ) — then the solution u(t,x) exists on the time interval [Xb,0) 
and can be uniquely extended to [Tq,0] as an element of (7°([To, 0], H k ~ 1 ) satisfying 


and 


||«||i°°(['Zb,0],VlX 1 .°°) < R , 


HP^mW- 

for To < t < 0. 


\\Pu(t)\\ 


P ± u 


H k-1 


H k ~ 1 


if k- (84/71 > 1 
if k /3 4 /7i = 1 > 
ifk- ( 84/71 < 1 

ifk- (84/71 > 1 or [B 0 ,] 


—t + (— t) 2 ( K ^ 4 / 11 ) ifk — fti/ji<l 



= 0 


Proof. First, the existence, uniqueness and continuation statements follow from standard results; for 
example, see [141 Ch.16 §1]. Therefore given a solution u = u(t,x) defined on the time interval [To,T„), 
T» £ (To, 1), to (1B.1MB.2I) . we act on (IB.II) on the left by D a B~ 1 to obtain 

B°d t D a u + B I d I D a u = jBD a Pu - B[D a , B^B^dtu - B[D a ,B~ 1 B I ]d I u + BD a (B~ 1 H). 

Using ED, we can write this as 


B°d t D a u + B I d I D a u= j BD a Pu - B[D a ^B^B^BPu^ + B[D a , B~ L B^B 0 )- 1 B 1 d lU 

3-inOnoOi-i 


1 B°](B°)~ 1 B I c 

1 lj 


-B[D a ,B- 1 B 1 }d I u + BD a {B- l H). (B.ll) 


As far as energy estimates are concerned, the only potential problematic term is the one with the coefficient 
j that becomes singular in the limit t 0. Since the rest of the terms can be estimated using well known 
techniques, again see |14[ Ch.16 §1], we will only estimate this potentially singular term in any detail. 

Setting a = 0 in (IB. 111) . we obtain, in the usual fashion for symmetric hyperbolic systems, the energy 
estimate 

—d t (u\B°u) = ~{u\BPu) + —(u\doBu) + — {u\ div Bu) + (u\H) + (u|G) (B.12) 

where 

div B = diB 1 . 

Defining the energy norm 

Nils = E ( D a u\B°D<*u ), 

|o;|<s 

we obtain from (IB. 12f) and the bounds (|B.3|) and (|B.10|) the estimate 
?k 

NHIIo < T l||PM|||N(«l^ 0 M}+7i||divi3||^|||M|||N2^(ll^llL 2 + ||G|| L Olhlllo^ T 0 <t < T*. (B.13) 
Using (IB.Ill , we can write 

(u|<9t-B°u) = (u\dtB°u) + -(m| [ D u B°(t , •, v, u ) • ( B°(t , ■, v, u))~ 1 B(t, ■, v , m)Pm]m) 


(B.14) 
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where 

d t B° = (, x , v, u) + D v B°(t, x, v, u ) • d t .v 


+ D u B°(t, x , v, u) 


(B°(t, x, v, u)) 1 


B 1 (t, x , v, u)diu + H(t, x , v, u) + G(t, x)) 


Observing that 

(u\ [D U B° ■ (B°)- 1 BFu]u) = ( F ± u\F- l [D u B° ■ (S°)- 1 SPu]P j -P- l m) + (P- L u\P ± [D u B° ■ ( J B°)- 1 BPu]PPu) 

+ ( Fu\P[D u B° ■ (B°)- 1 SPu]P- l P- l u) + < Fu\F[D u B° ■ (B°)“ 1 BPm]PPu), 

we see from (IB.5MB.8ll that 

\(u\ [D U B° • (B°)- 1 BPu]u>| < 4|t|0||u||£ a + P\\Fu\\ 2 L2 , (B.15) 

where /3 = J2i=i Pi- From (IB. 131) . (IB.14I) and (IB.lSj) . we see with the help of (1B.3I) that 

&IMIS < 2(ft -^ /7l) |p u |g+ 7l (461+11^511^ + 11 di y g|| £ -)| U Hg+2V ; n(||g||L2 + ||G||La)Hu|o (B. 16 ) 

for T 0 <f<T*. 

Next, using [P, B] = 0, we observe, for |a| < k, that 

(D a u\B[D a ,B~ 1 B 0 ](B 0 )~ 1 BFu) = (D a P- L u\B[D a ,B- 1 P ± B 0 P}(B 0 )- 1 BPu) 

+ (D a F- L u\B[D o ‘,B- 1 F ± B 0 ]F ± (B 0 )- 1 FBFu} + (D a Pu\B[D a ,B- 1 B 0 ](B 0 )- 1 BPu). 

Applying the Sobolev, commutator, product and Moser calculus inequalities, i.e. Theorems IA.21 [A.31 and 
IA.4I to the above expression, we find with the help of (IB. 41) and (IB. 91) the estimate 

(D a u\B[D a , B~ 1 B°](B°)~ 1 BPu) < C(K 1 ,K 2 )\\Fu\\ H 4Fu\\ H k-i 

where 

-fro = IMIl“([T 0 ,T,),.FT) and -fr 2 = H^Hl 00 ([To,0),// fc ) • 

Hence, for any e > 0, we find that 

(D a u\B[D a , B~ 1 B°](B°)~ 1 BFu) < C(K U K 2 ) (e||Pu||^ + c(e)||Pu||| 2 ) (B.17) 

by Ehrling’s lemma, Lemma IA.5I and Young’s inequality in the form ab < -y<z 2 + | b 2 for a, b > 0 and 
r > 0. We also see from 


(D a u\ [D U B° ■ (B 0 )~ 1 BFu\D a u) = 


:s,± D a u\P ± [D u B° ■ (B°)- 1 HPm]P- l P- l L> q u) + (F- l D°‘u\F ± [D U B° ■ (B°)- 1 BPu\PPD a u) 

+ (PD a u\P[D u B° ■ (B^BPulF^D 01 ^ + (PD a u\P[D u B° ■ (B°)- 1 BPu\PPD a u), 

the bounds (IB.5MB.8I) . and the calculus inequalities that 

\(D a u\[D u B°- (B°)- 1 BPu]l> a u)| < A\t\9\\D a u\\l 2 

+ C’(A' 1 ,cc)(/3i||Pu||| 2 + (ft + p 3 )\\Fu\\ L 2 \\FD a u\\ L 2 ) + 2/3 4 ||P-D“ U ||| 2 . 

Applying Young’s inequality to the this expression, we see for any e > 0 that 

\(D a u\[D u B° ■ (B^BFu^u)] < A\t\6\\u\\ 2 Hk +C(K 1 ,e^,co)(e\\Fuf Hk + c(e)||Pu||£ a ) +2f3\\PD a u\\ 2 L 2. 

(B.18) 

Setting 

v = l|G'llL°°([T o ,0),H fc ) 

and applying the standard energy estimates, i.e. as in the derivation of (IB.12I) . for symmetric hyperbolic 
systems to (IB. Ill) , we obtain, after summing over \a\ < k and using (IB.171) and (IB. 181) together with the 
calculus inequalities, the estimate^ 


d t M\l < '' IWull + CiKuK^K^B^P^v) 


for To < t < T*, where 


(e|P«lfc+c(e)||P«lo) +H 


k + v 


(B.19) 


K 3 — ||9t?l||L=o([ Toi o) 



Here, we are taking essentially the same approach to deriving energy estimates for symmetric hyperbolic systems as in 
Ch.16 §1]. 
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Since k — B Ik i > 0 by assumption, we see, after choosing e small enough and adding a suitable multiple 
of dBTfil) to dHUU), that 

9t (| U || + c|u|g + v 2 - J t ^||Pu(r)|| dr^j < C{K 1 ,K 2 , K 3 ,0, 71 ,/3, w, i/)(||u|| \ + v 2 ) 

where c = c(Ki,K 2 , K 3 ,9,k, w, v). Applying Gronwall’s inequality yields 

\\u{t)\\ 2 Hk + u 2 - I -||Pu(r)||^ fc <Ce c{t ~ To) (\\u(T 0 )\\ 2 Hk +v 2 ), T 0 < t < T*, 

Jt 0 t 

for some constant C(Ki,K 2 , K 3 , 9 ,71,72, fi, k, w, v). From this, Sobolev’s inequality and the continuation 
principle, it follows, for any fixed R > 0, that there exists a 

S = < 5 (i?, K- 2i K 3 , 9 , 71,72, /?, k, uj) > 0 

such that the solution u = u(t, x) exists on the time interval [To, 0) and satisfies 

rO 


\ u \\L°°([T o ,0),H k ) + v + 


1, 


- / -||ru(r;n H k 

JT 0 t 


2 dr ) < R 


(B.20) 


provided that ||u(To)|| ff fc + v < 5. 

To complete the proof, we assume that the condition ||zi(T 0 )|| ff fc + v < 5 holds. Writing EH as 

d t u = (B 0 )- 1 (B 1 di'u + -BPu + H + G 


and multiplying on the left by P 1 -, we obtain 

dtP^u = P ± {B°y 1 ( B : diu + H + G) + ^P ± (B°)~ 1 PBPu. 
Integrating this in time, we see after applying the H k ~ 1 norm that 
-P- L u(fi)|| i?fc -i < [ ||P- L (B°(r))" 1 T / (r)c) 7 u(r)|| fffc -i 

Jt 1 


(B.21) 


1, 


+ ||P- L (B u (r))- i U(r)||^-i + ||P- L (3? u (r))- i G(r)||^-i - -||P^( J B u (r))- i POT M (r)||^-i dr 

(B.22) 

for any t\,t 2 satisfying T 0 < t\ < t 2 . Using the calculus inequalities, i.e. product, Sobolev, Holder, and 
Moser, in conjunction with (IB.91) and the energy estimates (IB.201) . we obtain from the above inequality 
that 

\\P ± u{t 2 )-P ± u(t 1 )\\ H k-i < C(R, I\ 2 ) (|t 2 —U|—^ ^||Pu(t)|| Hk-idr^j < C{R, A" 2 )(|t 2 -ti|+o(|t 2 -U|))- 

It follows directly from this inequality that lim t/ *o P ± w.(t) exists in H k ~ l { T"), and moreover, that 

P- L U e C°([T 0 ,0],JT fc_1 ). 

Next, we apply the H k ~ 1 norm to (IB.211) and then integrate in time. With the help of the calculus 
inequalities, (IB. 91) . and (IB.20I) . this yields the estimate 


[ \\d t P ± u(T)\\ H k-i dr < C(R,I\ 2 ) T 0 < t < 0. 
Jt 0 


Multiplying m on the left by P shows that Pu satisfies 


B^Pu = jBPu + H + G 


(B.23) 


(B.24) 


where 


B l = PB l l 


B = P£>P, G = PG 


and 


H = PH + PB^d.P^u. 
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Applying the same arguments used to derive the energy estimate (IB.ldll to the equation (IB.241) yields 

| Pm I o 

(B.25) 


^|||Pm||| 2 < K ^ 4/7l |P«|g + ^(0+ \\d t BO\\ L «, + ||divB|U-)|Pu|S ++ IIGlMMo 


t 


for To < t < 0 , where now the energy norm ||| • ||| is defined by 

INIs = E ( D a u\B°D a u }, 


| q;|<s 


and 


d t B° = (<9 t B°)(t, x , v, u) + D v B°(t , x, v, u) ■ d t v 

+ D u B°(t, x, v, u) ■ ( B°(t , x , v, m)) _1 ( — B 1 (t, x, v, u)diu + H(t, x, v, u) + G(t, x)) 

Using (IB.4I) and (IB.20I) . we see that (IB.25I) implies that 

NMo < * ~ ^ hx 1 Pm 1 o + C(R,K 2 ,K 3 ,0, 7 i,^)[(l + H^mH^-i)|Pm || 0 + \\u\\ m + ||G || L2 
for T 0 < t < 0. It then follows from Gronwall’s inequalitjf] and the bounds (IB.201) and (IB.231) that 

—* if k — , 84/71 > 1 

tln(A-) if k — ^ 4/71 = 1 , To < t < 0. (B.26) 

4/Tl if /s - ^ 4/71 < 1 

Proceeding as above, we can apply the same arguments used to derive the energy estimate (IB.191) to 
(IB.241) to obtain the estimate 

1 


■u(t 


< 


^IPMlti < * P t 4hl |Pu| ti + C(R, K 2 , K 3 ,6 , 7l , v) 


Choosing e small enough, we see that the inequality 
d t \\\ru 


' ( e ll|P M lllfc - 1 + c ( e )l|P‘ u llo) + IP^IIU-i 

(B.27) 


1 , 


llk-i < f 4/7l) |||PM||k-i + C(R, K 2 , K 3 , e, 71,72, V, 0 4, k) ( l\\Pu\\ L 2 + l 


t 

holds for all a € (0,1), which, in turn, allows us to conclude via Gronwall’s inequality by choosing a close 
enough to 1 and using the estimate (IB.26|) that 






if k - , 84/71 > 1 


Tq < t < 0. 


(_f)<dK-/3 4 /7i) if a - / j 4 / 7 1 < 1 ’ 

From this estimate and (IB. 271) . we obtain after another application of Gronwall’s inequality that 

(-t if k - ,84/71 > 1 

||PM(t)|| ff fe—1 < < tln(^) ific —/ 3 4 / 7 ! = l , T 0 < t < 0 . 

71 if a- ^4/71 <1 

Using the above estimate in conjunction with (IB.22I) . we see with the help of the calculus inequalities 
and the bound (IB.201) that 

|*2 - *i| if k - 84/71 > 1 


m(* 2 )-P M(*i)|| ff fc-1 < 


1*2 - *i| + (-* 1 ) 2(R -^/^ l} - (-f 2 ) 2 ( s -^/^i) if ft - 84/71 < 1 


(B.28) 


for any *i,f 2 satisfying To < *1 < * 2 < 0. We also note in (IB.221) that if [T°,P] = 0, then the term with 
the - vanishes and we obtain the estimate 


i(* 2 ) — P^m(*i)|| j/fc —1 < C(R , A 2 )|* 2 — *i|, To < *1 < *2 < 0. 


(B.29) 


"'Here, we are using the following form of Gronwall’s inequality: if x(t) satisfies x'(t) < a(t)x(t) + hit,), t > To, 
then x(t ) < x(To)e A ^ + e A ^~ A G^h(r) d,T where A(t) = J 7 / a(r) dr. In particular, we observe from this that if 
a(t) = j + b(t), A G K, where \b(t)\ < r then 

x(t) < e^~ T °) x( To) (^) A + e r d- 2T °) ( —t) A £ dr 


for Tq < t < 0. 
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Sending t 2 Z' 0 in (IB. 281 ) and (IB. 291 ) . we see that P^m satisfies 


i(t) - 


H k ~ 1 


< 


-f 




if R — Pul'll > 1 or [T°,P] = 0 
if ii- (84/71 < 1 


for Tq < t < 0. 


□ 


Remark B.2. 

(i) The bounds ([B.3MB.8I) and (IB. 101) do not need to hold for all (v, u) £ K M x R N . The conclusions 
of Theorem IB.II remain valid provided that 

(a) there exists constants r, 1Z > 0 such that (IB.31) - (IB.81) and (IB. 101) hold for all (t,x,v,u) £ 
[To, 0] xT”x B r (R M ) x B-jz(R n ), 

(b) the map v £ C°([T 0 , 0], H k (T n )) satisfies |M|t~([T 0 ,o]xR") < r, 

(c) and the continuation principle is modified to having the solution u(t, x) remain in a proper 
subset of Bk(M. n ) for (t, x ) £ [To,T*) x T" in addition to satisfying |M|l°°([To,t„),w 1 ’“>) < 00 , 

(d) and the constant R > 0 is chosen so that R < TZ. 

(ii) Given any two constants 81)82 > 0, we observe that we can satisfy 

„ „ 28i , ^ ^ 282 

Gl < -r—: (77 and 62 < — , = 

w + IP-Lu ) 2 v / wT]P t m | 2 

for all u £ Biz(R n ) by setting 


81 = 


pi{n + n 2 ) 


81 = 


82 V(K + K 2 ) 


and u> = 1Z. 


This shows that if the bounds 


[ D u B°(t , x,v,u) ■ ( B°(t , x,v,u)) 1 B(t,x, v, M)Pw]P -L | op < \t\Q + 8 i|Pw| 2 , (B. 30 ) 

(P^ [D u B°(t, x,v,u) ■ ( B°(t, x,v,u))~ 1 B(t,x, v, u)Pm]P| op < \t\6 + 82|Pw| ; (B. 31 ) 


and 


\¥[D u B°(t, x, v, u) ■ (B°(t, x, v, u)) 1 B(t,x,v,u)Pu]P J ~\ o < \t\9 + 83 |Pw| 


(B. 32 ) 


hold for all ( t,x,v,u ) £ [T o , 0 ] x T” x B r (R M ) x Bk(M. n ), then the bounds (IB. 51 ) - (IB. 61 ) will be 
satisfied for all (t,x,v,u) £ [T o , 0 ] x T” x B r (R M ) x Bn(R N ) for the constants Pi,p2,P3 and u 
given by 


81 


Pi(n + K 2 ) 

2 


83 


82 


2 


and ui = 7 Z. 


In particular, this shows that by choosing 1Z small enough, we can always arrange that 81 > 82 
and 83 are as small as we like. 

Since the left-hand side of (IB. 6 I) - (IB. 8 I) are linear in Pm, it is clear that there always exists 
constants 82,83)84 > 0 such that the bounds (IB.31I) - (IB.32D and (IB. 81) hold for all (t,x,v,u) £ 
[To, 0 ] x T” x B t (R m ) x Bk(R n ). Moreover, by choosing 1Z > 0 small enough, we can take 84 > 0 
as small as we like. 


Thus we can conclude the following: 

If there exists a constant 81 > 0 such that 

|P^ [D u B°(t, x, v, u) ■ ( B°(t , x, v,u))~ 1 B(t , x, v, m)Pu]P _ l | oj) < |t|0 + 8 i|Pm | 2 

for all (t,x,v,u) £ [To,0] xPx B r (WL M ) x Bn(R N ) and 7 Z £ (0,7^o], then there exists constants 
8 i) 82 , 83 ) 84 > 0 such that the bounds (IB. 51) - (IB. 81) hold for all (t, x, v, u) £ [To, 0]xPx T r (R M ) x 
Bji(R n ). Moreover, by choosing 7 Z > 0 small enough, the constants 81)82,83)84 can be taken 
arbitrarily small. 
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